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Abstract 

In a full QCD lattice study with A'^^ = 2 Wilson fermions, we seek to optimize 
the signals for the disconnected contributions to the matrix element of flavour- 
singlet operators between nucleon states, which are indicative for sea quark effects. 
We demonstrate, in form of a fluctuation analysis to the noisy estimator technique, 
that - in order to achieve a tolerable signal to noise-ratio in full QCD - it is 
advantageous to work with a Z2-noise source rather than to rely only on gauge 
invariance to cancel non-gauge-invariant background. In the case of the ttN cr-term, 
we find that 10 Z2-noise sources suffice on our sample ( about 150 independent 
QCD configurations at /? = 5.6 on 16^ X 32 with 

^sea = 0.157, equivalent to 
M^^/Mp = 0.76(1)), to achieve decent signals and adequate fluctuations, rather 
than 300 such sources as recently used in quenched simulations. 



1 Introduction 



The direct computation of full fledged hadronic matrix elements containing flavour- 
singlet operators A = qTq, such as the vr-N cx-term, cttv, or the singlet axial vector 
current forward matrix element between nucleon states, has posed serious problems to 
lattice gauge theory. As a result there is up to now only faint evidence for sea quark 
effects from full QCD lattice computations of such matrix elements|l|]. 

The practical bottleneck is given by the computation of disconnected diagrams, i.e. 
ubiquitous insertions of A into quark loops disconnected from the valence quark graph, 
to say the nucleon propagator, P, in the vacuum background field. Technically, the 
disconnected insertion T = J2x[TrTM~}^^] must be calculated configurationwise, i.e. in 
correlation with the hadronic propagator. But the cost to compute all quark loops M~}^^ 
individually grows with volume which renders a direct evaluation of J2x prohibitively 
expensive 0. 

In the framework of quenched QCD, two groups |, 0, have recently tackled the 
problem by applying two different variants of the noisy estimator technique proposed 
some time ago for computing the chiral condensate X0- The trick is to induce the 
extensive character of the quark-loop insertion by inverting the Dirac operator M on a 
source extending over the entire (spatial) volume. This estimator technique is biassed 
in form of non-closed, i.e. bilocal and therefore non-gauge-invariant contributions that 
need be controlled. 

The Kentucky group ||, D pioneered the use of a stochastic volume source with Z2- 
noise (SETZ) to get rid of the 'bilocal crossterms', x y. On their small sample of 50 
configurations at /? = 6.0, they used about 300 such noise sources per configuration, at 
the expense of a substantial overhead in their computatinal cost. 

The authors of Ref. 0, on the other hand, utilize a homogeneous volume source 
(VST) and appeal to the gauge fluctuations to suppress the non-loop bias contributions 
via the Elitzur theorem. This approach however adds a substantial amount of noise to 
the data and may therefore require a very high statistics ensemble average to yield a 
reasonable signal. 

In exploratory quenched applications, where chromofield configurations can be cre- 
ated at low cost, both approaches have led to rather encouraging results on signal 
quality [0]. But in the context of a full QCD sea-quark computation, computer time pre- 
vents us at present from doing very high statistics sampling of gauge fields. So the ques- 
tion is whether effects from disconnected diagrams are accessible at all in small/medium 
size sampling of full QCD. 

This motivates us to look more deeply into the systematics of the stochastic estimator 
technique (SET) on actual full QCD configurations with Wilson fermions. To be spe- 
cific, we will proceed by studying the quality (size and fluctuations) of the signal from 
disconnected contributions to the scalar density. The strategy in applying estimator 
techniques should be such that - with minimal overhead - both systematic bias (from 
remaining nonlocal pollution) and fluctuations of the signal should become negligible 
compared to the statistical accuracy attainable from the size of the field sample. 

On the basis of the current statistics from our hybrid Monte Carlo productionp], ^ 
we will find that it is realistic to exploit the stochastic estimator technique for the study 
of quark loop effects in full QCD. 



2 Asymptotic Expectations 



In order to determine the disconnected part of one needs to calculate the expectation 
value (P Tr(M~^)), where P denotes the (momentum zero) proton correlation function 
and M is the fermion matrix. 

With VST one calculates ^i^^iC) on each gauge configuration C by solving 

M{C)X = h , (1) 

where h is a. volume source vector with components hi = 1. The average over the gauge 
configurations 

I Nc 



A=^j:nC)Y.Mr,HC) (2) 

^^c c=l i,j 

can be separated into local and non local contributions 



-t Nc T Nc 

^ = 1^ E PiC) E M-l{C) + — E PiC) E Mr.i(C) . (3) 

^^C c=l i c=l i^j 

As the latter are not gauge invariant they vanish (only) in the large Nc limit: 

lim A = {PTr{M-^)) . (4) 

Nc^oo 

In the stochastic estimator technique, one uses (complex) random source vectors with 
the property 

1 Ne 

lim ^T.v:{E,C)v,iE,C) = 5., (5) 

Ne-*oo Me E=i 

and computes the quantity rj^E, C)M~^{C)rj{E, C) repeatedly (A^^; times) on each con- 
figuration. 

It appears natural to choose rj{E, C) according to a Gaussian distribution (SETG) 
0. A decomposition similar to eq.(]^) then yields 
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where 
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Unfortunately the term Ton{NE,C) on the right hand side of eq.(^ introduces a gauge 
invariant bias, whose suppression definitely requires the number of estimates per config- 
uration to be in the asymptotic regime, where eq.^ becomes valid. 

This dangerous bias is removed from the beginning, if one samples the components 
of the random vector ri{E, C) according to a Z2 distribution. In this case, the relation 
T]* {E , C)rii{E , C) = 1 holds for any single estimate E, and Ton vanishes on each gauge 
configuration. 

For the remaining nonlocal bias, Toff{NE, C), SETZ combines two additive suppres- 
sion mechanisms: (a) the gauge fluctuations cancel them as non gauge invariant objects 
on a sufficiently large sample of gauge fields, even for A^^; = 1; (b) the Z2-noise kills 
them in the large Ne limit, even on a single gauge configuration. 

In order to understand the efficiency of the competing methods for rendering good 
signals in more detail, we consider the variance of {PTr{M~^)) in each case. Asymp- 
totically one finds 
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o'gauge ^he Variance calculated with the exact value of TrM~^ on each configuration, 
and a'^ff and a^^ are the variances due to the distribution of To// and To„ within the 
process of stochastic estimation. The cornered brackets stand for the average over gauge 
configurations. The abbreviation GOV in eq.|3 reads in its full length 
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for VST 



These formulae]^ show how the signals will be blurred by the fiuctuations of the estimates 
Tr{M~^{C)). Note that the terms in the large curly brackets of eq.^ are all gauge invari- 
ant. Hence they survive even in the limit of an infinite sample of gauge configurations 
unless suppressed otherwise : there is a way to fight them within SETZ and SETG ac- 
cording to a suppression^, but there is no way to infiuence them at all within VST. 



So far the discussion is qualitative only, as we do not know the relative size of these 
additional terms in eq.^ Furthermore, the calculation of requires the investigation of 
the ratio (P(0 t)Tr{M-^)) / {P{0 t)) rather than (P(0 t)Tr\M-^)) itself. This 

similar analysis has been carried out in for a modified version of VST, applied to the 
calculation of the tttt scattering length. 

^The terms depicted by COV behave like ^JI/Ne for SETZ and SETG. In case of SETG COV 
contains a gauge invariant part due to To„, which can be removed only in the limit Ne — > 00. All other 
terms in COV vanish for Nc — > 00, as they are not gauge invariant. 



can entail cancellations due to additional correlations between numerator and denomi- 
nator. Last not least we have to keep in mind that these considerations are only valid 
in the asymptotic limits of the gauge and estimator samples. In the next section we will 
therefore present a numerical study of the situation, under actual working conditions of 
full QCD. 

3 Numerical Results 

Our present analysis of af}''^'^ is based on 157 configurations from our ongoing Hybrid 
Monte Carlo run0, Q with two degenerate flavours of dynamical Wilson fermions. Here 
we work with Ksea = -157, which amounts to a ratio M^^/Mp = 0.76(1), and to a quark 
mass of rriq ~ LSm^. 

1. To set the stage, we consider first the implications of the stochastic estimator 
technique on the chiral condensate, x = Tr{M~^) /12V and compare the performance of 
Gaussian and Z2-noises in full QCD. In Fig.^ we show - on a given configuration - the 
standard error resulting from the two methods as a function of the number of estimates, 
Ne- Clearly SETZ is superior (by about a factor two) in obtaining a good signal for 
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Figure 1: The standard error of x for SETG and SETZ on a given configuration. 

X; we note that this is similar to the related quenched situation^]. For this reason we 
pursue SETZ in the following. 




2. An obvious way to economize is to properly adjust the accuracy in solving III 



eq.0. We therefore compute in the next step the chiral condensate in its dependence on 
the inversion accuracy r = || MX — h || / || X ||. Fig.|^ shows the convergence behaviour, 
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Figure 2: The difference Ax(r) = xi'^^max = 10 ^) —x{^) on a given configuration. The 
dashed horizontal lines indicate the one a margin from SETZ with Ne = 300. 

Ax = x(10^^) ~ The two horizontal lines in the plot refer to the la margin from 

the stochastic noise with Ne = 300 estimators, on a single field configuration. It can be 
seen that below r ~ 10~^, Ax is safely inside this margin. It is therefore sufficient to 
operate with the convergence condition r ~ 10~^. 

3. The next question relates to the size of the sample of gauge configurations required 
to observe a signal of the scalar density matrix element {P\qq\P). This quantity^ is 
extracted from 

(P(0 ^ t)Tr{M-^)) 



disc 



{Tr{M-')) const + t {P\qq\P)Tc 



(10) 



(P(0 - t)) 

To start we choose a large number of stochastic sources, Ne = 300 as proposed for 
quenched simulations in Ref. @], in order to avoid bias from the nonlocal pollutions to 
the trace. In Fig.^ we show how the linear rise in 

jK^fyisc evolves 

more and more clearly 

as the sampling is increased from 50 to 100 and 157 gauge configurations. The linear 

•^Throughout this work we appUed gauge invariant gaussian smearing |12 with N — 50, a — 4 to the 
Proton operator at the source. 
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Figure 3: i?(t)'^*"= measured with SETZ (Ne = SOO) on 50, 100 and 157 gauge field 
configurations. Linear fits were performed in the range t = 7 to t = 15. 
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Figure 4: measured with VST on 50, 100 and 157 gauge configurations. Linear 

fits were performed in the range t = 7 to t = 15. 



fits to this data yield the values 5.89(2.37), 3.44(1.45) and 2.51(0.77) for the slope. 

We conclude that it appears mandatory to work on samples of at least 100 configura- 
tions. Under this condition we (a) retrieve a reasonable signal and (b) find the statistical 



errors on the data points to follow the expected behaviour 1 / ^/Nl. For comparison we 
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Figure 5: TrM ^, R'^'^^'^{t = 11) and the corresponding variances as a function o/I/Ne- 
The data is taken from the full set of 157 gauge configurations. 

show in Fig.H R{tY^^'^ calculated with VST on 50, 100 and 157 gauge configurations. 
Note that the signal to noise ratio is much worse in this case. For the slope we find 
7.13(4.01), 5.40(2.8) and 3.50(2.20), consistent with SETZ, although with much larger 
statistical errors. 

4. Given the signal on our present sample of 157 configurations we can now ask 
the question whether and how far we can actually relax in the number of stochastic 
estimates, A^^;, without deteriorating its quality. The starting point is the observation 
that the fiuctuations from SET add to the inherent fiuctuations from the gauge field 
sample (see eq.H). For economics Ne should be chosen just large enough to suppress 
this undesired effect. Fig.^ shows the response of TrM~^ and of i?'^**'^(t = 11) as well as 
the variances a'^{TrM~^) and a'^{R'^^^'^{t = 11)) with respect to changes in 1/Ne- 
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Figure 6: R{t)'^'^^'^ measured on 157 gauge configurations, with 1, 10 and 300 stochastic 
Z2 estimates per configuration. 



We find that the mean values are not affected by the variation of Ne within the 
statistical accuracy. The variances however display a significant decrease over 1 < Ne < 
10, which ends in a somewhat noisy plateau. Obviously it does not pay to increase A'^^ 
beyond a value of 10. For comparison we quote also the results obtained with VST on 
the identical sample of gauge configurations: TrM~^ = 451817, R'^^'"^{t = 11) = 43.07, 
a'^{TrM-^) = 121863(14582), and a^{R'^''%t = 11)) = 75265(15073). Based on the 
variance we find SETZ to outperform VST substantially. In terms of statistics we 
recover a gain of a factor 2.5 to 3 from using SETZ instead of VST, in the present 
application. The full t-dependence of R'^^'^^{t) is displayed in Fig.^, showing again that 
10 estimates within SETZ are sufficient to produce a reasonable signal on our sample 
of 157 configurations, while the signal is nil with one estimate only. On the other hand 
300 estimates are definitely (and fortunately!) unnecessary. 



4 Conclusion and Outlook 

We have presented a full QCD study on the stochastic estimator technique applied to 
the disconnected diagrams of the scalar density. We found that it is indeed possible 
to achieve clean signals from only (9(10) Z2-noise sources with a weak convergence 
requirement, r ~ 10~^ on the iterative solver to TrM~^. This saves a factor 30 in 
computer time compared to previous applications of this technique in the quenched 
applications. 

With as few as 157 dynamical field configurations we obtain a reliable signal on 
< P\qq\P >disc.- This opens the door to perform a detailed analysis of the tt-N a term 



as well as the axial vector matrix elements in the proper setting, i.e. without recourse 
to the quenched approximation. Work along this line is in progress. 

It is obvious that in full QCD simulations the use of appropriately optimized Z2- 
noise techniques will be of utmost importance when it comes to the estimate of more 
complex quantitites in flavour- singlet channels, where the underlying correlators contain 
two disconnected fermion loops. 

Exploratory quenched applications of VST to compute the ?7'-mass and the annihi- 
lation diagrams in vrvr-scattering in the isospin zero channel at threshold have attained 



good signals [13, 14]. This makes us hope that medium size sampling will lead in full 



QCD to reliable signals as well, once optimized estimator techniques will be applied to 
grasp the loops. 
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